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Upper limits on the flux at the Sudbury Neutrino Observatory have been set based on the 
j/e charged-current reaction on deuterium. The reaction produces a positron and two neutrons in 
coincidence. This distinctive signature allows a search with very low background for i>e's from the 
Sun and other potential sources. Both differential and integral limits on the i/e flux have been 
placed in the energy range from 4 - 14.8 MeV. For an energy-independent — > conversion 
mechanism, the integral limit on the flux of solar Pe's in the energy range from 4 - 14.8 MeV 
is found to be < 3.4 x lO^cm'^s"^ (90% C.L.), which corresponds to 0.81% of the standard 
solar model *B flux of 5.05 x lO^cm^^s"^, and is consistent with the more sensitive limit from 
KamLAND in the 8.3 - 14.8 MeV range of 3.7 x 10^cm"^s"^ (90% C.L.). In the energy range from 
4-8 MeV, a search for Ve's is conducted using coincidences in which only the two neutrons are 
detected. Assuming a spectrum for the neutron induced fission of naturally occurring elements, 
a flux limit of <I>i,^ < 2.0 x 10^cm"^s"^ (90% C.L.) is obtained. 

PACS numbers: 26.65+t,13.15+g,14.60.St,13.35.Hb 
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I. INTRODUCTION 

This paper presents results from a search for i/g's 
with the Sudbury Neutrino Observatory (SNO) via the 
charged-current reaction (CC) on deuterons : 

De + d^e^ + n + n- 4.03 MeV. 

The distinctive signature of a positron in coincidence 
with two neutrons aUows SNO to search for D^s with 
very low background. By means of (n, n)-coincidence de- 
tections, SNO has sensitivity to i>e's with energies above 
the reaction threshold of 4.03 MeV. For coincidences in- 
volving a e+, SNO is sensitive to Pe's above a threshold 
of 4 MeV -I- E*|,"'^;[, where ^l^laii is the analysis threshold 
applied to the recoil positron. We present results for in- 
tegral De flux limits in the neutrino energy range from 4 - 
14.8 MeV under the assumption that j>e's originate from 
a ®B spectrum, and below 8 MeV under the assumption 
that DgS originate from a fission spectrum. Differential 
limits on the De flux have been placed, independent of 
any particular spectral assumptions. 

As a solution to the solar neutrino problem [J, anal- 
yses of global solar neutrino data favor matter-enhanced 
neutrino oscillations with mixingparameters in the Large 
Mixing Angle (LMA) region 0,13 ■ Studies of solar neu- 
trino data have demonstrated that approximately two- 
thirds of ®B solar neutrinos convert to active flavors other 
than before reaching Earth HHIllllSll. The deficit 
of reactor antineutrinos reported by the Kamioka Liquid 
scintillator Anti-Neutrino Detector (KamLAND) experi- 
ment Jl supports the LMA solution to the solar neutrino 
problem under the explicit assumption of CPT conser- 
vation in the neutrino sector. More generally the spin 



flavor precession (SFP) mechanism H US 

or neutrino 

decay could contribute to the observed neutrino fla- 
vor transformation at a sub-dominant level by converting 
a small fraction of solar I'e to Dg,. In SFP models, neutri- 
nos are assumed to have a transition magnetic moment 
of the order of 10^^^ Bohr magnetons [l3 if they are of 
Majorana type. Solar magnetic fields, which are known 
to be time dependent, couple to this magnetic moment to 
convert v,, into with a combination of neutrino flavor 
oscillations and SFP mechanisms. Neutrino decay mod- 
els allow heavier neutrino mass eigenstates to decay into 
light D^. 

Nuclear fission Vf. spectra peak at low energies and fall 
approximately exponentially with energy, and have neg- 
ligible intensity above 8 MeV J^]. Positrons produced 
in CC reactions from fission v^s are too low in energy 
to be detected by the present analysis. However, by 
conducting an analysis involving only (n, n)-coincidences 
it is possible to study the energy region from 4-8 MeV, 
providing some sensitivity to v^s that might originate 
from naturally occurring neutron induced fission sources. 
Because the fission spectrum is significant in this energy 
region, whereas the *B spectrum is peaked at higher ener- 
gies, a separate analysis of this region for j/g's originating 
from fission is performed. The expected flux from man- 
made reactors would provide a negligible contribution to 
this analysis. 

Additional sources of D,, 's are cosmic ray interactions in 
the upper atmosphere and the diffuse supernovae back- 
ground. The flux of these types of v^s are dominated 
by Ve energies above the 4 - 14.8 MeV energy range 
investigated in the present analysis and hence are only 
addressed peripherally by the differential limits below 
15 MeV. 
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II. EXPERIMENTAL DATA 

SNO is an imaging water Cherenkov detector located 
at a depth of 6010 m of water equivalent in the Inco, Ltd. 
Creighton mine near Sudbury, Ontario, Canada. SNO 
detects neutrinos using an ultra-pure heavy water target 
contained in a transparent acrylic spherical shell 12 m in 
diameter. Cherenkov photons generated in the heavy wa- 
ter are detected by 9456 photomultiplier tubes (PMTs) 
mounted on a stainless steel geodesic sphere 17.8 m in 
diameter. The geodesic sphere is immersed in ultra-pure 
light water to provide shielding from radioactivity in both 
the PMT array and the cavity rock. The SNO detector 
has been described in detail in [T3|. 

The data reported here were recorded between Novem- 
ber 2, 1999 and May 28, 2001 and span the entire first 
phase of the experiment, in which only D2O was present 
in the sensitive volume. In comparison to earlier SNO 
analyses 0, 01 j a cut employed to remove events follow- 
ing muon candidates was modified to maximize rejection 
of false Df, candidates. Any event with an assigned ki- 
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netic energy above 18 MeV and all events following that 
event within 0.5 s are rejected. These modifications re- 
sulted in a total livetime of 305.9 live days, a reduction of 
less than 0.2% from SNO's neutral-current (NC) analysis 

H 

III. ANALYSIS 

Interactions of v^. 's with deuterons produce a positron 
and two neutrons. The positron generates a prompt 
Cherenkov signal, while the neutrons must first thcrmal- 
ize before generating 6.25 MeV gamma rays from their 
capture on deuterons. The mean neutron capture time is 
42 ms and the diffusion length is ~ 110 cm in pure D2O. 
Inside the D2O volume Ve events can be identified by a 
coincidence of 2 or 3 particles. 

The analysis procedure consists of two steps. The first 
is similar to the SNO solar neutrino analysis and is de- 
scribed in 0. In this step, PMT times and hit pat- 
terns are used to reconstruct event vertices and assign 
to each recoil positron or electron a most probable ki- 
netic energy, Tcff. The recoil threshold in this analysis 
was Teff > 5 MeV, providing sensitivity to positrons and 
neutrons from the CC reaction. No i^e's originating from 
the conversion of solar i/g's are expected to have energies 
in excess of 14.8 MeV. Hence, the 18 MeV upper limit 
on recoil positron or electron energy does not remove po- 
tential solar or fission De candidates. A fiducial volume 
was defined by requiring reconstructed event vertices to 
be within 550 cm of the center of the acrylic vessel. This 
reduces both the number of externally produced back- 
ground events, and the systematic uncertainties associ- 
ated with optics and event reconstruction near the acrylic 
vessel. 

The second step of the analysis identifies coincidences 
among the accepted events. The size of the coincidence 
window, chosen to be 150 ms, was optimized to maximize 
the sensitivity to 2-fold coincidences in the presence of 
the background of accidental coincidences. Antineutrino 
candidates that are part of a burst of 4 or more events 
of any energy are discarded. Even under the assumption 
that 100% of the solar are converted to j/g's, four- 
or higher- fold coincidences are 10'^ times more likely to 
originate from a background such as atmospheric i^'s or 
spontaneous fission of ^38^ than from a 3-fold event 
in coincidence with a random event, as determined by 
Monte Carlo simulations. 



A. Detection Efficiencies 



ration, and correctly describes the detector status as a 
function of time. Based on the number of simulated and 
extracted coincidences inside the signal region, the 2- and 
3-fold coincidence event detection efhciencies were found 
to be: 

e(e+,„,n)= 1.11 ±0.02(stat.);H2(syst.)% 
e(e+,n) = 10.27 ± 0.05(stat.) to:94(syst.) % 
e(„,„) = 1.20 ± 0.02(stat.) ^HoCsyst-) % 

Among coincidences originating from Vf, interactions, 
(e"*" , n)-coincidences are 10 times more likely to be de- 
tected than (n, n)-coincidences or (e"'", n, n)-coincidences. 
While the (n, n)-coincidences can originate from v^s 
with energies as low as 4 MeV, coincidences contain- 
ing a positron must originate from D^s with energies 
above 9 MeV. Figure ^ displays the number of expected 
Vf, events as function of the antineutrino energy under 
the assumption that i>eS originate from a *B spectrum 
with a total flux of 5.05x 10^cm~^s~"'^. The above coinci- 
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Neutrino Energy (MeV) 

FIG. 1: Number of expected coincidence events as a func- 
tion of the antineutrino energy, for Vs's originating from a *B 
spectrum with a total flux of 5.05x 10^cm~^s~^. Coincidences 
containing a positron and (n, n)-coincidences are shown as 
dashed and dotted histograms, respectively. The solid line 
represents all types of detected coincidences as a function of 
Ve energy. 



Coincidence detection efficiencies were determined 
from a simulated Vf, data set. The Monte Carlo (MC) 
simulations sampled VeS from a *B energy spectrum 
[1^ with a total flux 100 times the standard solar 
model (SSM-BPOO) flux [H. The simulated data 
set matches the experimentally recorded data set in du- 



dence detection efficiencies are consistent with the aver- 
age detection efficiency estimated from Monte Carlo 
ee+ = 40.09 t\f2 (syst) ± 0.10 (stat)%, and an average 
neutron detection efficiency of e„ = 14.38 ± 0.53% jg. 
The neutron response and systematic uncertainty on the 
response were calibrated with a ^^^Cf source. 
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B. Backgrounds 

Backgrounds to the search can be divided into 
two categories: coincidences caused by S'e's from known 
sources, and coincidences from processes other than Vf, 
charged-current reactions. These are presented in Ta- 
ble n The first category has contributions from atmo- 
spheric, reactor, and diffuse supernovae i^e's. The back- 
ground contribution from atmospheric v^^ is estimated 
to be 0.07 ± 0.01 coincidences. It is derived from the 
CC cross sections coincidence detection efficiencies, 
and a parameterized Vf^ spectrum extrapolated into the 
energy range below 50 MeV [13. Highly energetic at- 
mospheric i^e's only contribute to the background if de- 
tected as (n, n)-coincidences, due to the applied upper 
bound on the recoil lepton energy. The background es- 
timate from nuclear power reactors yields 0.019 ± 0.002 
coincidences. The calculation is based on the monthly 
reported actual power output of all commercial reactors 
[so] within 500 km of SNO and an average reactor v,, 
spectrum |12L . Furthermore, the Ve flux is assumed 
to be reduced as a result of neutrino oscillations gov- 
erned by the best fit oscillation parameters tan^(6')=0.41 
and Am2=7. 1x10-5 eV^ 0. The number of CC in- 
teractions from diffuse supernovae neutrinos was esti- 
mated by combining a theoretical flux calculation 
with SNO's detection efficiency and contributes < 0.005 
coincidences. The calculation is consistent with experi- 
mental limits jl^l . Electron antineutrinos from the decay 
chain of terrestrial radioisotopes such as ^38^ and ^^^Th, 
geo-antineutrinos, do not contribute to the background 
because their energies are below the threshold of the CC 
reaction. 

The main contributions to the second category of back- 
grounds originate from atmospheric possible ^ssy gg_ 
sion events in the detector media, and accidental coin- 
cidences. The atmospheric v background is estimated 
to account for 1.46^0 45 coincidences. Atmospheric neu- 
trinos were sampled from a realistic spectrum |l7j and 
propagated with the neutrino-interaction generator NU- 
ANCE '2l'| . Detailed event information, including energy 
and multiplicities at the depth of SNO, was then pro- 
cessed by a full detector simulation. 

Spontaneous fission background from ^38^ detec- 
tor media was estimated to contribute less than 0.79 co- 
incidences. This limit was derived from an inductively 
coupled plasma mass spectrometer (ICP-MS) measure- 
ment of the concentration of 66±17 (fg U)/(g D2O) 
[8ll |. for which approximately 14 spontaneous fission de- 
cays of 238]j are expected per kton-year. The detector 
response was calculated on the basis of detection effi- 
ciencies and a discrete probability distribution for neu- 
tron multiplicities in spontaneous fission 22] . The above 
estimate is an upper limit since the 238u concentration 
measurement was obtained on September 4, 2003, after 
the data set reported here, and after the addition of NaCl 
to the heavy water. Because the reverse osmosis purifi- 
cation system, which was operated while the data accu- 



mulated for the present analysis, could not be run with 
NaCl in the heavy water, the measurement is likely to 
include more ^38^ contamination than was present when 
the data were taken. 

Accidental coincidences are formed by individual 
events that pass the analysis cuts and have a random 
time correlation. Their number has been calculated as a 
function of the time-dependent singles rate in the detec- 
tor and amounts to 0.13 ^0 94 coincidences. This number 
was confirmed independently by a direct measurement of 
signals following within 150 ms of a large number of ran- 
dom triggers. 

The background from neutron capture reactions on 
oxygen, which can produce multiple gamma-rays above 
the deuteron photodisintegration threshold, was esti- 
mated to be less than 0.05 coincidences. It was calculated 
from abimdances of ^^O and ^*0, relative intensities of 
gamma rays produced in ^^0(n,7)^^0 and ^^0(n,7)^^0 
reactions, the ratio Pi7o,iso/-Pd of the neutron capture 
probability on ^''O and ^^0 and deuterium, as well as 
the total number of observed neutrons. Instrumental 
backgrounds are events produced by electrical pickup or 
emission of non-Cherenkov light from detector compo- 
nents. Their background contribution is determined to 
be < 0.027 coincidences (95% C.L.). The number of coin- 
cidences of this type is assessed by means of a bifurcated 
analysis, which employs sets of orthogonal cuts aimed 
at instrumental background rejection. The background 
from a-capture reactions on carbon, which can produce 
a neutron in coincidence with an — e~ pair, was found 
to be less than 1.7x10^^^ (90% C.L.) coincidences. It was 
estimated on the basis of the total number of neutrons 
in the signal region and a MC calculation. Other back- 
grounds originate from radioisotope contamination and 
can produce coincidences through /3-7 or /3-rt decays, but 
are found to be entirely negligible. They are estimated 
on the basis of their respective radioisotope contamina- 
tion levels. 

Upper limits and uncertainties on individual back- 
grounds have been combined under the assumption that 
they are independent. The total background amounts to 
1.68 to;45 coincidences. 



IV. RESULTS 

The search for v^. candidates in the experimental data 
set employs the same cuts on energy and fiducial volume 
as well as the same coincidence extraction algorithms as 
were used to derive the Monte Carlo-based coincidence 
detection efficiencies. One 3-fold and one 2-fold coinci- 
dence were found. Table ITU summarizes the characteris- 
tics of the two candidate coincidences and their con- 
stituent events. On an event-by-event basis it is not pos- 
sible to uniquely identify individual constituent events as 
a positron or a neutron. Therefore, the analysis re-groups 
(e^ ,n)- and (n, n)-coincidences into a single category of 
2-fold coincidences. This category has an order of mag- 
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TABLE I: Types of coincidence backgrounds and number 
of expected coincidences in the SNO detector for the data 
set. Upper limits and uncertainties on individual backgrounds 
have been combined under the assumption that they are in- 
dependent. 



Type of z-'e 



background 



expected 
coincidences 



Atmospheric 
Reactor 

Diffuse supernovae 
geo-antineutrinos 



0.07 ± 0.01 
0.019 ± 0.002 
< 0.005 
0.0 



Total ; 



background 



0.09 ±0.01 



Non-, 



background 



expected 
coincidences 



Atmospheric u 


1 46 

-0.45 


^^^U spontaneous fission in detector media 


< 0.79 


Accidental coincidences 


n 1 +0.06 
^■^^ -0.04 


''0(n, 7)^+10 , wiiere x=17,18 


< 0.05 


Instrumental contamination (95% C.L.) 


< 0.027 


^"Cia.ne+e-y^O (90% C.L.) 


< 1.7 xlO"^ 


Intrinsic: 




^^^Bi: /3-7 decay 


7.6x10"^ 


^^°T1: 0-n decay 


« 10"** 


^"•^Tl: /3-7 decay 


8.7x10"'' 


7 — » Compton -f photo-disintegration n 


< 8 xlO"" 


Total non-Pe background 


1 59 +"-'^^ 


Total background 


1 fis 

l-o8 45 



TABLE II; Two z/g candidate coincidences are found. Listed 
are kinetic recoil lepton energy and radial position for each 
constituent event as well as spatial separation and time sep- 
aration relative to the first particle in each coincidence. 





candidate 


Toff (McV) 


'"fit (cm) 


Ar (cm) 


At (ms) 


I 


1st particle 


8.58 


283.2 


0.0 


0.0 




2nd particle 


5.39 


472.4 


206.7 


16.7 




3rd particle 


5.15 


349.2 


178.3 


20.3 


II 


1st particle 


6.95 


506.4 


0.0 


0.0 




2nd particle 


6.09 


429.5 


81.8 


88.9 



nitude higher sensitivity than 3-fold coincidences and is 
sensitive to t'e's with energies as low as 4 MeV. 



A. Differential Limits 

This analysis sets model-independent differential lim- 
its on the Df, flux in the neutrino energy range from 4 
- 14.8 MeV. Bin sizes in neutrino energy were chosen to 
be 1 MeV. Based on the observed 2-fold coincidence, and 
under the conservative assumption of zero background, 
the Bayesian upper limit on the number of 2-fold 
coincidences amounts to 3.89 at the 90% C.L. In each 
neutrino energy bin, it is assumed that the candidate 2- 
fold event was produced by an De of that energy. The 
upper limit on the number of candidate events is then 
corrected for detector acceptance and cross section to ob- 
tain a limit on the absolute Pe flux at that energy. As a 
result, the limit in each energy bin is model-independent 
and maximally conservative, but limits in different en- 
ergy bins are strongly correlated. Systematic uncertain- 



ties in the theoretical cross sections, energy resolution 
differences between data and MC, simulation failures, as 
well as systematics related to data reduction combine 
to about 2% and have been taken into account. Figure 
13 displays Pe flux limits for the energy range from 4 - 
15 MeV at greater than 90% C.L. Super-Kamiokande's 
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FIG. 2: Limits on the />e flux from SNO (bars) and SK (stars). 
Bars represent 90% C.L. flux limits for 1 MeV wide energy 
bins and are based on the assumption that the observed 2-fold 
coincidence originates from that particular energy bin. Stars 
indicate limits for mono-energetic u^'s . 

(SK) flux limits for mono-energetic De's are shown for 
comparison [241] . Super-Kamiokande's limits are based on 
data events, after subtraction of spallation background, 
which fall in the ±lcr range of a Gaussian that describes 
the detector response to mono-energetic f/g's. The SNO 
and SK limits are slightly different in nature since SNO 
limits were calculated for a series of 1 MeV wide bins in 
neutrino energy. 



B. Integral Limit 

Under the assumption that the energy distribution of 
solar Pe's follows a ^B spectrum, and that both observed 
candidates are of solar origin, an integral limit on the 
solar De flux is derived. The 2- and 3-fold coincidences 
are joined in order to maximize the sensitivity. Using 
an extended Feldman-Cousins method .25, 26] to include 
the background uncertainty in the form of a two-sided 
Gaussian, the 90% C.L. upper Hmit for 2 candidate co- 
incidences and 1.68 ILq'^ background events has been 
calculated to be 3.8 coincidences. A MC calculation was 
used to convert a given flux into a number of ob- 
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served events. The 3.8 coincidences translate into a 
flux limit of $ < 3.4 x lO^cm^^s""'^ in the energy range 
from 4 - 14.8 MeV. The systematic uncertainties have 
been treated similarly to the differential analysis. The 
analysis energy window contains 83.4% of the SSM-BPOO 
Ve flux of 5.05lJ:^i X lO^cm-^s-i [l^. The above 
limit implies a 90% C.L. upper bound on the conversion 
probability of solar i/e's into Pg's of 0.81%, if are 
assumed to follow a ^B spectrum. This assumption is 
equivalent to an energy-independent conversion 
mechanism. 

If the analysis is restricted to the energy range from 
4-8 MeV, only the observed 2-fold coincidence repre- 
sents a Ve candidate since the 3-fold coincidence could 
only have originated from a Ve with an energy in excess 
of 12.6 MeV. Within the 4 - 8 MeV energy window the 
background is conservatively assumed to be zero coinci- 
dences. Using a Bayesian prescription [23, the 90% C.L. 
upper limit on one candidate and zero assumed back- 
ground corresponds to 3.89 events. Assuming a fission 
spectral shape 0, 0| from possible naturally occurring 
elements, this defines a limit of $ < 2.0x 10^cm~^s^^ in 
the energy range from 4-8 MeV which covers 9% of the 
average reactor v,, flux. 



C. Comparison with other Experiments 

Previously, other experiments have set very stringent 
limits on the flux from the Sun. Under the assump- 
tion of an unoscillated *B spectral shape for solar v^s 
KamLAND [23| limits the solar flux of Ve& to less than 
3.7 X 10^ cm-^s"^ (90% C.L.). This measurement is 
based on the neutrino energy range from 8.3 - 14.8 MeV 
and corresponds to an upper limit on the Vf, — > 
conversion probability of 2.8 x 10"'' (90% C.L.). 

SK's integral flux limit is based on the energy range 
from 8-20 MeV and, under the assumption of a ^B 
spectrum for solar i?e's and a total flux of 5.05 x 
10^ cm~^s^^, places a 90% C.L. limit on the conversion 
probability of less than 0.8% (2^. This corresponds to 
an absolute flux limit of 1.4 x 10^ cm~^s~^ in the energy 
region from 8-20 MeV which contains 34.4% of the 
total ®B flux. The present SNO analysis investigates the 
energy range from 4 - 14.8 MeV and uses a completely 
independent direct counting method with very low 
background. Therefore it can set a comparable limit on 
the Ve flux despite the fact that SK's exposure is 800 
times larger. Due to the time dependence of the solar 
magnetic field the v^. flux originating from conversion 
of solar ^B neutrinos could vary as function of time. 
Table Hill specifies the existing limits on the flux and 
indicates the time frames during which the relevant data 
were recorded. 

The present analysis provides a 90% C.L. upper limit 
for Ve energies from 4-8 MeV of 2.0x10^ cm^^s^^ 
assuming a neutron induced flssion spectral shape. 
Because the contribution from man-made reactors in the 



TABLE III: Integral flux limits !?e/SSM Ve at the 90% 
C.L. relative to the SSM-BPOO *B flux jlj, periods of time 
during which the relevant data were taken, and energy ranges 
over which the various experiments have sensitivity to solar 
Ve are presented. 



Experiment 


Time Period 


Energy 


(MoV) Limit(%) 


KamLAND I27| 


3/4/2002 - 


12/1/2002 


8.3 


- 14.8 


0.028 


SNO 


11/2/1999 


- 5/28/2001 


4.0 


- 14.8 


0.81 


SK [241 


5/31/1996 


- 7/15/2001 


8.0 


- 20.0 


0.8 


LSD [281 


before 4/1996 


8.8 


- 18.8 


1.95 


Kamiokandc [291 


6/1988 - 4/1990 


12.0 


- 13.0 


5.07 



region of the SNO detector is calculated to be very low, 
(see Table nj, this can be considered as an upper limit on 
the flux from naturally occurring neutron induced flssion 
sources. However, we note that if all 54 events observed 
by KamLAND for Ve energies from 3.4-8 MeV were 
considered to originate from the neutron induced fission 
of naturally-occurring elements rather than from nearby 
man-made reactors, an upper limit on the flux in 
the 3.4 - 8 MeV energy region of 1.55 x 10^ cm~^s~^ 
(90% C.L.) can be derived for the KamLAND location. 
During the completion of this paper, KamLAND 
reported a new result of 258 events in a 515-day live 
time and a 33% larger fiducial volume. No significant 
change in the Ve flux limit is expected as a result. 

The neutron detection efficiency in pure D2O was 
14.4% 0. For the phase of the SNO experiment in which 
NaCl was added to the D2O, improving the neutron 
detection efficiency to 39.9% |8|, an increase in the ef- 
fective Ve sensitivity by a factor of about 3 is anticipated. 



D. Neutrino Decay 

Although flavor transformation of solar neutrinos is be- 
lieved to be dominated by matter-enhanced neutrino os- 
cillations with mixing parameters in the LMA region, 
SFP mechanisms or neutrino decay may also contribute. 
Since neutrino decay is expected to be energy-dependent, 
SNO's low Ve energy threshold of 4 MeV is a valuable fea- 
ture to test for non-radiative neutrino decay of the form 
V2 — > vi + X. Here V2 and vi refer to the heavier and 
lighter neutrino mass eigenstates and X is a scalar parti- 
cle (e.g. a Majoron) IjjJ . For quasi-degenerate neutrino 
masses, a lower limit on the lifetime T2 of the heavier 
neutrino is found to be T2/m2 >0.004 s/eV. For hierar- 
chical neutrino masses, the limit amounts to T2/m2 > 4.4 
xlO"^ s/eV, equivalent to T2 > 0.44 ^s if TO2 « 0.01 eV 
[3^. Previously, KamLAND |27| has presented lower lim- 
its on non-radiative neutrino decay based on the Ve en- 
ergy range from 8.3 - 14.8 MeV, and found T2/TO2 >0.067 
s/eV for quasi-degenerate and r2 > 11 /is for hierarchical 
neutrino masses. 
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V. CONCLUSION 

In summary, our analysis represents a novel detection 
technique to search for C'e's, with very low backgrounds. 
Based on the one 2-fold and one 3-fold observed coinci- 
dence, integral limits on the i^e flux in the energy range 
below 8 MeV and in the range from 4 - 14.8 MeV have 
been set under the assumption of a fission and a *B spec- 
trum, respectively. Spectrally independent differential 
limits have been placed as well. The derived limit on the 
flux of solar Pe's was used to constrain the neutrino life- 
time. Within SNO's sensitivity we independently confirm 
the SK and KamLAND results on fluxes. 
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